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We extend the e-expansion of continuum chiral perturbation theory to nonzero lattice spacing in 
the framework of Wilson Chiral Perturbation Theory. We distinguish various regimes by defining 
the relative power counting of the quark mass m and the lattice spacing a. We observe that 
for m ^ aAQf^Q, the explicit breaking of chiral symmetry in Wilson fermions is still driven by the 
quark mass and lattice corrections are highly suppressed. The lattice spacing effects become more 
pronounced for smaller quark masses and may lead to non-trivial corrections of the continuum 
results at next-to-leading order We compute these corrections for standard current and density 
correlation functions. A fit to lattice data shows that these corrections are small, as expected. 
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1. Introduction 

Matching results of light quarks lattice QCD simulations with the predictions of the chiral 
effective theory is an important step towards a comprehensive picture of low energy properties of 
strong interactions. Through this matching one can eventually extract the Low-Energy Couplings 
(LECs) of the effective theory: many results have been presented in the past months for the leading 
order (LO) and next-to-leading order (NLO) couplings, both for the SU(2) and the SU(3) effective 
theory (for a recent review see for instance [1]). While lattice QCD simulations in the chiral regime 
experience constant progress, the control over the systematic uncertainties remains a relevant issue. 
Lattice artifacts and finite-volume effects must be thoroughly understood. Moreover, since the 
matching is done at a given order in the chiral effective theory, the effect of the neglected higher 
order contributions must be under control. 

For this reason, it is very useful to extract the LECs from a large set of independent observables 
and from different kinematic regimes. An interesting approach is to investigate QCD in a finite 
volume V = Lr'T in the so-called £-regime [2], where the pion wavelength is larger than the size of 
the box, MjiL < L In this situation volume effects ai^e enhanced, while mass-effects are suppressed 
with respect to the usual infinite-volume case (or p-regime, where M-jiL » 1). The consequence 
is that, at a given order in the perturbative expansion, less unknown LECs will appear, making the 
£-regime potentially convenient for the extraction of the LO constants. 

Lattice simulations in the £-regime are challenging, since one still needs a fairly large volume 
(L 3> 1/ {AnF)), while the quark mass must be very small (mZV ^1)- This fact strongly influences 
the choice of the lattice action to be used. Ginsparg-Wilson (GW) fermions [3] have many theoret- 
ical advantages, since they guarantee exact chiral symmetry at finite lattice spacing [4]. The price 
to pay to maintain continuum-like chiral symmetry is the high computational cost. While a large 
number of quenched computations in the e-regime with Ginsparg-Wilson Dirac operator has been 
performed (see [5] for a recent study, and references therein for precedent computations), dynami- 
cal simulations are still very costly. For recent calculations with Nf = 2 dynamical chiral fermions 
see [6, 7, 8, 9]. H. Fukaya presented at this conference new preliminary results for A'y = 2 + 1 
obtained by JLQCD/TWQCD [10]. 

On the other hand, dynamical simulations with the Wilson Dirac operator with 0{a) improve- 
ment are becoming fairly inexpensive. Nevertheless, since chiral symmetry is explicitly broken at 
finite lattice spacing, the spectral gap is not bounded from below: at small quark masses the prob- 
ability distribution of the gap can approach zero values, leading to possible integration instabilities 
and sampling inefficiencies in the algorithm. In [11] the empirical stability bound m > m,„i„, with 
a/ \fV has been established. A possible solution to this problem comes from a reweighting 
technique [12] (see also [13]). This method has been applied in [14] to simulate Wilson fermions 
in the e-regime. Lattice simulations in the £-regime are feasible also with Wilson Twisted Mass 
fermions, as presented in [15] and also at this conference [16]. 

In principle, the matching with the chiral effective theory should be performed only after a contin- 
uum extrapolation of the lattice results. While this is not an unrealistic goal for the near future, the 
presently available simulations in the £-regime ai^e carried out at a single value of the lattice spac- 
ing. In [14] the pseudoscalar and axial congelations functions turned out to be very well described 
by continuum chiral effective theory at NLO. Still, it is important to have a theoretical understand- 
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ing of the impact of explicit breaking of chiral symmetry on computations in the £-regime. We 
address this question in [17]:the tool that we adopt is the so-called Wilson Chiral Perturbation The- 
ory (WChPT) [18, 19], the low-energy effective theory for lattice QCD with Wilson Dirac operator. 
A similar analysis has been canied out in [20, 16]. 



2. Wilson Chiral Perturbation Theory 

The chiral effective Lagrangian of WChPT is an expansion in powers of pion momenta p^, 
the quark mass m and lattice spacing a. Based on symmetries of the underlying Symanzik action 
[21], the chiral Lagrangian including all terms of 0{p'^,p^m,m^,p^a,ma) is given in [19]. The 
0{a^) contributions are constructed in [22, 23]. In the following we consider the case Nf = 2 with 
degenerate quark mass m. The leading order Euclidean chiral Lagrangian in the continuum is given 
by [24, 25] 

^2 = -jTr {d^Ud^U"^) - '-^Tv {U + f/t) . (2. 1) 

The pseudo Nambu-Goldstone modes are parametrized by the SU(2) field U {x) = exp {lit, [x]) /F, 
and F, B = L/F^ are the familiar LO couplings. The leading terms involving the lattice spacing 
are 

= dW45Tr (dfiUdfiU^) Tr {U + f/"^) - amWes (Tr {U + U"^))^, (2.2) 

^,2 = —C2a^Tr{U + U^)) , (2.3) 

where m = 2Bm and d = 2Woa. W45, Wes, Wq and C2 are new LECs. Note that the mass parameter m 
in Eq. (2.3) is the so-called shifted mass [18]: besides the dominant additive mass renormalization 
proportional to 1 /a it also contains the leading correction of 0{a). 

Currents and densities in WChPT can be constructed by a standard spurion analysis or by 
introducing source terms. Here we report the axial vector current and the pseudoscalar density 
including the leading 0{a) corrections [26, 27, 28]: 

^;;,wchPT = ^;;,cont { 1 + ^dW45Tr{U + [/^) I , P^chPT = ^ont { 1 + ^aWeHTr{U + 1/^ ) I , 

(2.4) 

where 

A;;,cont = yTr {T"{U^d^U - Ud^U')) , P,Vt = ''^Tr {T"{U - U^)) , (2.5) 
and T" are SU(2) generators normalized such that Tr{T"T^) = d"'' /2. 

2.1 Power counting in infinite volume 

In WChPT there are two parameters which break explicitly the chiral symmetry, the quark 
mass m (counted as 0{p^)) and the lattice spacing a. The power counting is determined by the 
relative size of these two pai^ameters. In particular, one defines [26, 29] two different regimes: (i) 
the GSM ^ regime, where m ~ aA^i^^ and (ii) the Aoki regime where m ^ a^AQ^.^. In the Aoki 

'GSM stands for generically small masses. 
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regime lattice artifacts are more pronounced, and the _Sf„2 in Eq. (2.3) enters already at LO. The 
pion mass at leading order is given by 

GSM regime : Ml = 2Bm, (2.6) 
Aoki regime : Mq = 2Bm - lc2a^ . (2.7) 

The sign of C2 governs the phase diagram of the theory [18]. 

2.2 Power counting in the £ -regime 

In the e-regime the chiral limit is approached by keeping jj. = niLV < 0(1); this corresponds 
to the situation where the pion wavelength is much larger than the linear size of the box. The 
main effect of formulating the effective theory in this regime is that the pion zero-mode becomes 
non-perturbative and its contribution has to be treated exactly. This is achieved by factorizing the 
pseudo Nambu-Goldstone boson fields as 

U{x) = exp(^j^{x)^Uo, (2.8) 

where the constant Uq G SU(2) represents the collective zero-mode. The non-zero modes ^ can 
still be treated perturbatively. The £-regime requires a reorganization of the perturbative series: in 
the continuum, this corresponds to taking the quark mass of order m ~ 0{e^). 

In order to extend the WChPT to the e-regime, one has to assign a relative power counting of 
the lattice spacing a with respect to the quark mass m. By assuming that the quark mass can be 
considered of order m ~ 0{e'^) also in WChPT ^, we obtain 

GSM regime : m ~ 0{aA^cD) ar^ 0[e^), (2.9) 
Aoki regime : m 0{a^ Aqq^)) — > a ~ 0{e^). (2.10) 

In addition, the e- expansion allows us to introduce an intermediate counting between the GSM 
and the Aoki regime: we can define the GSM* regime, where a 0{e^). 

We compute mesonic two-point functions within the WChPT in the £-regime. In particular, we 
give explicit results for the pseudoscalar and axial time correlators, 

5"'Cpp{t) = jd'x{P"{x)P\0)), d"'CAA{t) = j d'x{Al{^)4m- (2.11) 

Currents and densities are defined in Eq. (2.4); the subscript "WChPT" is now omitted. With our 
power counting, we find that in the GSM regime lattice coiTcctions enter only at NNLO, while in 
the Aoki regime effects of lattice artifacts show up already at LO. 

In the following we concentrate on the intermediate GSM* regime, where corrections appear 
at NLO. At this order, they are given only by the term in Eq. (2.3). This suppression comes 
from the fact that the lattice spacing corrections in the chiral effective theory action and in the 
effective operators are either quadratic in a or they come with an additional power of either m 
or p^. Notice that this is valid also for unimproved Wilson fermions: a non-perturbative 0{a) 
improvement removes corrections due to as well as the 0{a) terms in the operators, i.e. acts 
only on terms which are subleading in our power counting. 

^While this is a natural choice in the GSM regime, the situation in the Aoki regime can be more subtle. See [17] for 
a more detailed discussion on this issue. 
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2.3 Leading corrections in the GSM* regime 



The continuum pseudoscalar and axial correlators at NLO in the £-expansion can be written 
as [30] 

Cpp^AA ct(0 = ap^A + bpAh (t/T) , 



(2.12) 



with a pai^abolic time dependence given by the function hi (t) 
For Nf = 2 the coefficients ap^A,bpA explicitly read [30] 
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(2.14) 



V /i(2At) 

/i, /2 ai^e modified Bessel functions of the first kind; /3i and ^oo ai"e so-called shape factors [31, 30], 
which depend only on the geometry of the finite box. Egft is the quark condensate at one loop [2] 

3 Pi 



Teff = 11 + 



(2.15) 



2F2 

and /leff = mLsffV. As anticipated, the continuum NLO expressions contain only the LO LECs £ 
and F. The first non-trivial modification of the continuum NLO results appear in the GSM* regime 
and it is due to Jf^2 only. In this case we can write down the full NLO coiTclators in WChPT as 

Cpp,AA{t) = Cpp^AAct{t) +CppAAcfi{t)- (2.16) 

By performing the explicit computation (see [17] for the full details) it turns out that the corrections 
CppAA (0 are time-independent, hence affect only the constant part of the correlators. In particular 
we obtain: 



Cppa2 



where 



2 " " ' T 

5At/f (2ai) - lO/i (2^1)72(2^1) - 3^11(2^) 
2^3/2(2^) 



(2.17) 
(2.18) 



and p = F^C2a^V is the dimensionless LEG which parametrizes the 0{a'^) correction. 
It is useful to compute the leading 0{a^) corrections to the PGAG quark mass: 



mpcAC = m 



1+P 



At' 



M/2(2m)/ 



(2.19) 



It is now possible to express the con^elators Cpp^AA{t) as a function of p. = mpcAC^V; the result is 



3v2 



Cpp{t) = Cppct{t) + -^pAa2 , Caa {t 



f2 . 
CAAa{t) + YP^a^^ 



where the continuum correlators are as in Eq. (2.12), but with the replacements pL 
Aeff = wpcAC^^eff^ and 

~ ^ 4/|2/3(2/l) -nfill{2fi)h{2fi)+2{3 - 2fi2)/i(2At)/|(2A) + 5A/|(2/1) 

2m\mhm 

Other correlation functions can be computed along the same line. 



(2.20) 

Meff 

(2.21) 
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Figure 1: Fit of the WChPT predictions to lattice data. All data points within the 
fit range of t/a E [6,18] for the four sea quark masses are included in the combined 
fit. The hopping parameter K = (0.128150,0.128125,0.1281,0.128050) corresponds to ampcAC = 
(0.0019(4), 0.0024(3), 0.0030(3), 0.0044(3)) respectively. The axial vector coiTelator is multiplied by a 
factor 50 for better visibility. 



3. Reanalysis of lattice data and conclusions 

These predictions from WChPT at NLO can be tested against lattice data generated in [14], 
where pseudoscalai" and axial con^elators have been computed on an ensemble with Nf = 2 flavors 
of dynamical improved NHYP Wilson fermions [32]. At a lattice spacing a ~ 0. 115 fm, two lattice 
extents are considered, Li = 16a ~ 1.84 fm and L2 = 24a ~ 2.8 fm. Quark masses approach the 
£-regime, with p. ~ 0.7 — 2.9 for the volume Vi = L\ and /i ~ 2.1 — 5.0 for the volume V2 = L\. In 
the GSM* regime, we have only the additional LEG C2 with respect to the continuum case. Notice 
that its value will depend on the particular- discretized action which is used. We simultaneously fit 
the two correlators for all available quark masses; for the volume V2, a fit in the range t € [6, 18] 
gives 

= 249(4) MeV, F = 88(3) MeV, C2 = 0.02(8) GeV^ (3.1) 

The data, along with the theoretical curves, are shown in Fig. 1. The errors from the renormal- 
ization factors Za, Z^^{ij. = 2 GeV) computed in [14] are not included in the uncertainties of the 
LEGs. Varying the time range of the fit does not give significant differences for the LEGs, as long 
as fmin/fl' > 4. Also discarding the heaviest mass does not change the results of Eq. (3.1) within 
the statistical errors. The smallest volume Vi yields values which are consistent with Eq. (3.1), but 
the large of the fit may indicate that NLO formulae are no longer applicable. The values of F 
and £ are compatible with other determinations [1], while the value of C2 is compatible with zero. 
A continuum fit (with C2 = 0) yields virtually unchanged values for F and Z, showing that cut-off 
effects do not impact the extraction of the LEGs beyond the level of the statistical uncertainties. 
This is a very encouraging result: simulations with Wilson fermions in the £-regime are feasible 
and seem to be a viable alternative to dynamical simulations with GW fermions. Lattice computa- 
tions on a wide range of lattice spacings and volumes would of course be very useful to test if the 



E^^(At =2GeV) 
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predicted NLO scaling is verified. The results derived here can be generaUzed in various ways, for 
example to the case with a twisted mass term or to an arbitrary number of flavors. 
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